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Abstract 

We show that spherical accretion onto astrophysical black holes can 
be considered as a natural example of analogue system. We provide, for 
the first time, an exact analytical scheme for calculating the analogue 
Hawking temperature and surface gravity for general relativistic accre- 
tion onto astrophysical black holes. Our calculation may bridge the gap 
between the theory of transonic astrophysical accretion and the theory of 
analogue Hawking radiation. We show that the domination of the ana- 
logue Hawking temperature over the actual Hawking temperature may be 
a real astrophysical phenomena. We also discuss the possibilities of the 
emergence of analogue white holes around astrophysical black holes. Our 
calculation is general enough to accommodate accreting black holes with 
any mass. 

1 Introduction 

By providing an analogy between the propagation of light in a curved geometry 
and the acoustic propagation in a steady, inviscid, transonic barotropic fluid, it 
has been demonstrated El E El El that the representative fluid posses the 
signature of an acoustic black hole (BH hereafter) embedded by an acoustic hori- 
zon at its transonic point (where the fluid Mach number becomes unity). This 
acoustic horizon emits thermal radiation known as analogue Hawking radiation. 
The temperature of such radiation is called analogue Hawking temperature, or 
acoustic Hawking temperature (Unruh 1 was the first to propose and compute 
such a temperature). Throughout our work, we use Tah to denote the analogue 
Hawking temperature, and Th to denote the the actual Hawking temperature 
where Th = hc^ /SnkhGMBH, Mbh being the mass of the black hole. We use 
the words 'analogue', 'acoustic' and 'sonic' synonymously in describing the hori- 
zons or black holes. The theory of analogue Hawking radiation opens up the 
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possibility to experimentally verify some basic features of BH physics by creat- 
ing the sonic horizons in the laboratory. A number of works have been carried 
out to formulate the condensed matter or optical analogs of event horizons [7]- 
It is also expected that analogue Hawking radiation may find important uses in 
the fields of inflationary models, quantum gravity and sub-Planckian models of 
string theory 

Since the publication of the seminal paper by Bondi in 1952 jS], the transonic 
behaviour of accreting fluid onto compact astrophysical objects has been exten- 
sively studied in the astrophysics community (including the fast ever paper on 
curved acoustic geometry by Moncrief and the pioneering work by Unruh 
in 1981 ^ initiated a substantial number of works in the theory of analogue 
Hawking radiation with diverse fields of application stated above. No attempt 
has yet been made to bridge these two categories of research, astrophysical 
BH accretion and the theory of analogue Hawking radiation, by providing a 
self-consistent study of analogue Hawking radiation for real astrophysical fluid 
flows, i.e., by establishing the fact that accreting black holes can be considered 
as a natural example of analogue system. Since both the theory of transonic 
astrophysical accretion and the theory of analogue Hawking radiation stem from 
almost exactly the same physics, the propagation of a transonic fluid with acous- 
tic disturbances embedded into it, it is important to study analogue Hawking 
radiation for transonic accretion onto astrophysical black holes and to compute 
Tah for such accretion. 

In this paper we show that a spherically accreting black hole naturally provides 
a concrete example of analogue system found in nature. We consider general rel- 
ativistic (GR) spherically symmetric accretion of inviscid polytropic fluid onto a 
Schwarzschild BH of mass Mbh- We first formulate and solve the conservation 
equations for such accretion and then demonstrate that such flows may become 
transonic depending on initial boundary conditions. We then, for the first time 
^ , completely analytically calculate the exact location of the acoustic horizon 
[TT] and calculate the quantities on the acoustic horizon which are required to 
compute the surface gravity and Tah for our model. We then analytically cal- 
culate the exact value of Tah in curved space time for all possible real physical 
spherical accretion solutions. 

It is important to note that the accreting astrophysical BHs are the only real 
physical candidates for which both the horizons, actual (electromagnetic) and 
analogue (sonic), may exist together simultaneously. Hence our application of 
analogue Hawking radiation to the theory of transonic astrophysical accretion 
may be useful to compare some properties of these two kind of horizons, by the- 
oretically studying and comparing the behaviour of the same flow close to these 

^Note that Unruh's original calculation of T/^u assumes the positional constancy of 
sound speed, which may not always be a real physical situation. This was modified to include 
the position dependent sound speed by Visser |1] for flat space and by Bilic |^ for curved 
space. However, none of these papers provides the exact calculation of rj, and acoustic horizon 
related quantities from first principles, and hence the numerical value of Tah in o-'^U existing 
literature (even for flat space, let alone curved space) is obtained using approximate, order of 
magnitude, calculations only. 
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horizons. As an illustrative example of such a procedure, we also calculate and 
compare Tah and Th for the same BH and for the same set of initial boundary 
conditions, and show that for astrophysically interesting regions of parameter 
space, Tah may well exceed Th- One should, however, note that both Th and 
Tah comes out to be quite less compared to the macroscopic classical fluid tem- 
perature of accreting matter. Hence BH accretion system may not be a good 
candidate to allow any observational test for detecting the analogue radiation. 
Our work rather introduces the most relevant classical system to manifest such 
analogues in curved space-time. Also we find the possibility that an actual BH 
event horizon may be embedded by an acoustic white hole. 



2 Formalism and results 



We take the Schwarzschild radius Vg = 2GMbh/c^, and scale the radial coordi- 
nate r in units of fg, any velocity in units of c, and all other physical quantities 
are scaled accordingly. We set G = c = Mbh = 1- The mass of the accreting 
fluid is assumed to be much less compared to Mbh (which is usually reality 
for astrophysical BH accretion), so that the gravity field is controlled essen- 
tially by Mbh only. Accretion is 6 and (j) symmetric and posses only radial 
inflow velocity. We concentrate on stationary solutions. Accretion is governed 
by the radial part of the GR time independent Euler and continuity equations 
in Schwarzschild metric. The conserved specific flow energy £ (the relativistic 
analogue of Bernoulli's constant) along each stream line reads £ = hut [T^ 
where h and are the specific enthalpy and the four velocity, which can be 
re-cast in terms of the radial three velocity u and the polytropic sound speed a 
to obtain: 
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7- (l + a2) \/ 1 ^.2- W 

In this work we concentrate on positive Bernoulli constant solutions. The mass 
accretion rate M may be obtained by integrating the continuity equation: 



M = Anpur 



I r-1 
r (1 — u"^) ' 



(2) 



where p is the proper mass density. A polytropic equation of state, p — Kp'' , 
is assumed (this is common in astrophysics to describe relativistic accretion) to 
define another constant S as: 



'E. — M — All pur 
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(7-1) 
7- {l + a?) 



(3) 



r(l-u2) 

Simultaneous solution of eq. (1-3) provides the dynamical three velocity gradient 
at any radial distance r: 

du _ u{l- u^) [a^ (4r - 3) - l] _ M 
~ 2r (r - 1) (u2 _ a'^) ^ V' 



(4) 
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The denominator of cq. (4) becomes zero at r = 1 (actual horizon) and when 
u = a (acoustic horizon). To maintain the physical smoothness of the flow 
for any r > 1, the numerator M and the denominator "D of eq. (4) have to 
simultaneously vanish. Hence by making A/" = I? = 0, wc get the so-called sonic 
point conditions (value of dynamical and sound velocities in terms of horizon 
location r/,) as 



Uh = ah = 



Substitution of Uh and ah into eq. (1) for r — provides: 



rl + r^Fi + rhT2 + Ts = 0, 



(5) 



(6) 



where 



Ti = 



2f 2 (2 _ 37) + 9 (7 - 1) 



4(7-1) (£2-1) 



To = 



£^ (37-2)^-27(7-1)^ 
32 (^2 - 1) (7 - if 



To = 



27 



64 (f 2 _ 1) • 



(7) 



Solution of eq. (6) provides the location of the acoustic horizon in terms of only 

two accretion parameters {5,7}, which is our two parameter input set to study 
the flow, and hereafter will be denoted by 7^2 ■ We solve eq. (6) completely an- 
alytically by employing the Cardano-Tartaglia-del Ferro technique. We define: 



griTs - 27r3 



2rf 
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e = cos-i 



Si 



Si 



Ef, n± = (i7i±i^2) 



(8) 



so that the three roots for Vh come out to be: 



-^-i(0_-^^/3n_). 



(9) 



However, note that not all V?i{i = 1,2,3} would be real for all V2- It is easy 
to show that if ^E" > 0, only one root is real; if '5 = 0, all roots are real and 
at least two of them are identical; and if \l/ < 0, all roots are real and distinct. 
Selection of the real physical [vh has to be greater than unity) roots requires 
the following discussion. 

Although we analytically calculate rh and other variables at r/j, there is no 
way that one can analytically calculate the flow variables at any arbitrary r. 
One needs to integrate eq. (6) numerically to obtain the transonic profile of 
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the flow for all range of r, starting from infinity and ending on to the actual 
event horizon. To do so, one must set the appropriate limits on V2 to model 
the realistic situations encountered in astrophysics. As £ is scaled in terms of 
the rest mass energy and includes the rest mass energy {£ = 1 corresponds to 
rest mass energy), £ < 1 corresponds to the negative energy accretion state 
where radiative extraction of rest mass energy from the fluid is required. For 
such extraction to be made possible, the accreting fluid has to posses viscosity 
or other dissipative mechanisms, which may violate the Lorenzian invariance. 
On the other hand, although almost any f > 1 is mathematically allowed, 
large values of £ represents flows starting from infinity with extremely high 
thermal energy, and £ > 2 accretion represents enormously hot (thermal energy 
greater than the rest mass energy) initial (r — *■ 00) fiow configurations, which 
are not properly conceivable in realistic astrophysical situations. Hence we set 
1 ^ £ ^ 2. Now, 7=1 corresponds to isothermal accretion where accreting 
fluid remains optically thin. This is the physical lower limit for 7; 7 < 1 is not 
realistic in accretion astrophysics. On the other hand, 7 > 2 is possible only for 
superdense matter with substantially large magnetic field (which requires the 
accreting material to be governed by GR magneto-hydrodynamic equations, 
dealing with which is beyond the scope of this paper) and direction dependent 
anisotropic pressure. We thus set 1 < 7 < 2 as well, so V2 has the boundaries 
1 ^ {£, 7} ^ 2. However, one should note that the most preferred values of 7 
for realistic BH accretion ranges from 4/3 to 5/3 \V6\. 

Coming back to our solution for r/i, we find that for our preferred range of 
7^2, one always obtains 'J < 0. Hence the roots are always real and three real 
unequal roots can be computed as: 




We find that for all I ^ V2 ^ 2, becomes negative. We then find that 
{V/j,^r/i} > 1 for most values of our astrophysically tuned V2- However, it 
is also found that \h does not allow steady physical flows to pass through it; 
either m, or a, or both, becomes superluminal before the flow reaches the ac- 
tual event horizon, or the Mach number profile shows intrinsic fiuctuations for 
r < rfi- We obtain this information by numerically integrating the complete 
flow proflle passing through Hence it turns out that we need to concentrate 
only on V?i for realistic astrophysical BH accretion. Both large £ and large 7 
enhance the thermal energy of the flow so that the accreting fluid acquires the 
large radial velocity to overcome a only in the close vicinity of the BH. Hence r^, 
anti-correlates with V2- To obtain (dv/dr) and {da/ dr) on the acoustic horizon, 
we apply L' Hospital's rule to eq. (4) and we have: 
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where 

$12 = -$2/2$l, $123 = ^$l-4$i$3/2$i, 
$1 = = [4r, (7 - 1) - (37 - 2)] 

$3 ^ ^(""^^1 [r^ (7 _ if - rn [lOl^ - 197 + 9) 
(4rh-3)5 

The GR description of the sonic geometry for a generahzed acoustic metric 
is provided by BiUc 0. Keeping in the mind that the surface gravity can be 
obtained by computing the norm of the gradient of the norm of the KiUing fields 
evaluated at r?i, we use the methodology developed in Bilic [HI to calculate the 
general relativistic expression for surface gravity, and Tah for our flow geometry 
is obtained as: 



he 

Tah 



inkbCMBH 



rl {rn ~ 0.75) 
{rh - 1)^ 



(13) 



where the values of rh, {du/dr)h and {da/dr)h are obtained using the system 
of units and scaling used in our paper. The stationarity approximation allows 
us not to discriminate between the acoustic apparent and event horizon. Note 
again, that, since and other quantities appearing in eq. (13) are analyti- 
cally calculated as a function of 1^2, we provide an exact analytical value of the 
general relativistic analogue Hawking temperature for all possible solutions of 
an spherically accreting astrophysical black hole system, something which has 
never been done in the hterature. If \/4r/i — 3(1/2— l/$4)($i2 — $123) > 1, one 
always obtains {da/dr < du/dr)h, which indicates the presence of the acoustic 
white holes at rh- This inequality holds good for certain astrophysically relevant 
range of V2] see following discussions. 

For a particular value of 7^2, we now define the quantity r to be the ratio of 
Tah and Th, i.e., t = Tah/Th- This ratio r comes out to be independent of 
the mass of the black hole, which enables us to compare the properties of two 
kind of horizons (actual and acoustic) for an spherically accreting BH with any 
mass, starting from the primordial holes to the super massive black holes at 
galactic centers. Note, however, that we applied the analogue to describe the 
classical perturbation of the fluid in terms of a field satisfying the wave equa- 
tion on a curved effective geometry. Our main motivation in this paper is not 
to rigorously demonstrate how the phonon field generated in this system could 
be made quantized. To accomplish that task, one can show that the effective 
action for the acoustic perturbation is equivalent to the field theoretical action 
in curved space, and the corresponding commutation and dispersion relations 
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Fig. 1: Variation of ratio of the analogue to the actual Hawking temperature 
(r) on conserved specific energy (£) and polytropic indices (7) of the flow. 

(see, e.g., may directly follow from there. However, such calculations are 
beyond the scope of this paper and would be presented elsewhere. 

Figure 1 shows the dependence of r on 7^2- The dependence is highly non- 
linear. We plot r corresponding to the black hole solutions only (solutions for 
which {da/dr)ii > [du/dr)^). Note that Tah > for some 7^2- In figure 2, 
we scan the complete V2 space by depicting its four important non-overlapping 
regions. For accretion with V2 taken from the lightly shaded regions marked 
by < Tah < Th, the Hawking temperature dominates over its analogue 
counterpart. For high 7, high £ flow ('hot' accretion), analogue Hawking radia- 
tion becomes the dominant process compared to the actual Hawking radiation 
{Tah > Th)- For low £ and intermediate/high 7, acoustic white holes appear 
(white region marked by WH, where da/dr < du/dr at r/j). For V2 belong- 
ing to this region, one obtains outflow (outgoing solutions with Mach number 
increasing with increase of r) only. This has also been verified by obtaining 
the complete flow profile by integrating eq. (6). The dark shaded region (lower 
right corner of the figure) represents V2 for which comes out to be physical 
{rh > 1) but $123 becomes negative, hence {du/dr)h and {da/dr)h are not real 
and Tah becomes imaginary. Note that both Tah > Th and white hole regions 
are obtained even for 4/3 < 7 < 5/3, which is the range of values of 7 for most 
realistic flows of matter around astrophysical black holes. Hence we have shown 
that the domination of the analogue Hawking temperature over the actual Hawk- 
ing temperature and the emergence of analogue white holes are real astrophysical 
phenomena. We did not perform formal stability analysis on white hole solu- 
tions. It may be possible that the white hole solutions are unstable, according 
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Fig. 2: Parameter space classification for P2 for four different sets of values of 

to the recent proposal that the sonic white holes posses intrinsic instability |15| . 
However, such detailed investigations are beyond the scope of this paper and 
we would like to defer it for our future work. 



3 Discussion 

In this paper we have studied spherically symmetric accretion onto astrophysical 
black holes as a natural example of analogue system. If accreting matter contains 
non-zero intrinsic angular momentum density, accretion becomes axisymmetric 
and forms disc like structure surrounding the black hole. Calculations presented 
in this paper may be extended for astrophysical BH accretion discs as well. Bilic 
in] studied the relativistic acoustic geometry for 2D axisymmetric space time and 
one can use such formalism to study analogue Hawking radiation for GR and 
post-Newtonian multi-transonic, rotating, advective BH accretion discs in the 
similar way as is done in this letter. Viscosity, however, is quite a subtle issue 
in studying the analogue effects for disc accretion. From the astrophysical point 
of view, even thirty years after the discovery of standard accretion disc theory, 
exact modeling of viscous transonic BH accretion, including proper heating and 
cooling mechanisms is still quite an arduous task, even for Newtonian flow, let 
alone for general relativistic accretion. And from analogue model point of view, 
viscosity is likely to destroy the Lorenzian invariance, and hence the assumptions 
behind building up an analogue model may not remain that much consistent. 
However, there are two possibilities through which the calculations presented 
in the paper may be translated onto the idea of disc accretion. Firstly, one 
can study the disc system in its weak viscosity limit. Secondly, if viscosity is 
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present, one can calculate the ratio of the in- fall to viscous time scale to see at 
which length scale the ratio comes out to be less than one. At that length scale, 
accretion acquires such a high radial in-fall velocity (measured on the equatorial 
plane of the disc) that the residence time of matter at a particular annulus of the 
disc becomes quite low and hence the angular momentum transport may become 
inefficient. So the flow may practically be taken as inviscid in that length scale. 
Now if it is found that the length scale at which the sonic horizon is formed, is 
more or less comparable or less than the length scale where viscosity becomes 
insignificant, then it should not be any problem in studying the analogue model 
even for disc accretion. The whole parameter space for which one can model 
the analogue process using accretion disc, will then depend on the accretion 
parameters for which the acoustic horizon forms at a distance equal to or less 
than the radial distance where the in-fall time scale is smaller compared to 
the viscous time scale. Our preliminary calculations show that indeed there 
are situations for which the above condition is self-consistently satisfied. Such 
calculations are in progress and will be presented elsewhere. 
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